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ABSTRACT 

A model for high temperature superconductors based on the idea of Cooper pairs 
comprised of electrons from different bands is studied. We propose that the two bands 
relevant for the cuprates are comprised of Cu d x 2_ y 2, d z 2, planar O p a , and apical O p z 
orbitals. Along the diagonal, k x = k y in the Brillouin zone, the two band Fermi surfaces 
may cross. We associate the optimal doping for the highest T c with this point because 
only in the vicinity of this touching point are inter-band Cooper pairs energetically possible. 
Due to the lack of time reversal invariance of an inter-band Cooper pair with itself, the 
standard interpretation of Josephson tunneling is altered such that the detailed nature of 
the single particle tunneling matrix elements contributes to the supercurrent. The d x 2_ y 2 
gap observations from Josephson tunneling are shown to arise from our model with pairing 
due to phonons. A Hubbard model is written down for the two bands at the Fermi energy 
with realistic parameters for Lai.ssSro.isCuC^. The anomalous normal state features in the 
nmr are calculated and qualitatively explained as due to the character of the two bands in the 
vicinity of the crossing point. The Hall effect is calculated using standard Bloch-Boltzmann 
transport theory. The observed strong temperature dependence of the Hall coefficient is 
reproduced and is due to the strong reshaping of the current carrying band Fermi surface 
due to band repulsion with the other band for dopings very close to the Fermi surface 
touching point. Reasonable quantitative agreement is also obtained for the nmr and Hall 
effect. A linear resistivity at optimal doping is expected due to the proximity of the second 
band in k space which can strongly relax the current and the "smallness" of the current 
carrying Fermi surface. 
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I. INTRODUCTION 

In a previous publication the author proposed an Inter-Band Pairing (IBP) model 
for high temperature superconductors. 1 We suggested that the fundamental idea of Cooper 
pairs composed of electrons in states k t and —k I be retained but rather than k t and 
—k I being from the same band, we consider pairs such that k t and —k I come from 
different bands at the Fermi surface. Such Cooper pairs, if they exist, are not time reversal 
invariant with themselves in contrast to all BCS-like (Intra-Band) Cooper pairs. The full 
Hamiltonian is time reversal invariant, though. We showed that with such Inter-Band pairs, 
the orbital character of the two bands plays a significant role in all Josephson tunneling 
experiments. The phase difference across a junction is no longer simply the phase difference 
of the two gap functions but also includes a contribution due to the phases of the hopping 
matrix elements across the junction. With BCS-like Cooper pairs, the product of the single 
particle tunneling matrix elements is mod-squared due to the time reversal invariance of the 
pair. This is not the case with IBP. We suggest that the observed 2 ' 3 d x 2_ y 2 character of the 
Josephson tunneling can be resolved by an "s-like gap" coupled with a u d x 2_ y 2" Cooper pair. 
These issues will be explained in detail below. 

Clearly, in order to have any possibility of creating an Inter-Band Cooper pair, 
{k ^(j, —k li) where U, L are labels for the two distinct bands (eventually taken to represent 

the Upper and Lower bands), the single particle energies ej^ = e^?, = must 
be close to the Fermi energy. Such a circumstance is easily available when the U and L 
band Fermi surfaces intersect. Of course, bands repel each other in general unless there is 
a symmetry that forbids mixing. Thus, Cooper pairing will occur in the vicinity of special 
symmetry points in the Brillouin zone. 

Our previous paper dealt with the unphysical situation of two distinct bands with 
exactly coincident Fermi surfaces at a given doping which was associated with the optimal 
doping for superconductivity. In this regard, the previous study of IBP must be regarded 
as an overly simplistic model that serves to illustrate some of the general principles of IBP 
but cannot account for quantitative features. This paper considers realistic bands for the 
cuprates and derives some quantitative consequences of the theory. 

We propose the two relevant bands for high T c are comprised predominantly of Cu 
d x 2_ y 2, Cu d z 2, planar O p a , and apical O p z orbitals. The two bands are strong mixtures of 
these orbitals except at special k points only. A simple Hubbard model is written down for 
the system. The parameters used in our Hubbard model are derived in the following paper 
by Perry and Tahir-Kheli 4 using ab-initio calculations on finite clusters for Lai .ssSro.lsCuO^ 
Contrary to the results of local density approximation (LDA) band structure calculations 
on the cuprates, we find two bands at the Fermi energy comprised primarily of the above 
orbitals. The difference between our results and previous band structure calculations are 
attributed to correlation effects. The optimal doping is the doping where the two band 
Fermi surfaces touch. This can happen along the diagonal where k x = k y . Laig5Sro.i5Cu04 
is primarily considered here so the 3D dispersion is small. 

The NMR spin relaxation rates at the planar Cu and O nuclei are estimated with these 
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bands and so are the respective Knight shifts. Reasonable quantitative agreement with many 
of the anomalous NMR properties 5-11 is obtained. The dramatically different temperature 
dependencies of the planar Cu and O spin relaxation rates is due to the conspiring of several 
important aspects of the orbital characters of the two bands U and L at the Fermi energy. 
The major reason is that there are two possible orbitals on the Cu, d x 2_ y 2 and d z 2, whereas 
there is only one important orbital p a on the planar O. The second most important reason 
is that the L band is almost full at optimal doping and the orbital character of the available 
L k states is dominated by the character at the top of the band at k — (it, it). 

The large spin relaxation anisotropy 10 on the Cu site is discussed and we argue this is 
due to a small mixture of d xy character in the two bands. We show that mixing of d xy of a 
few percent is sufficient to produce the observed anisotropy ratio of pz 3.4. 

The Cu and O Knight shifts are shown to have similar T (temperature) dependence to 
the O relaxation rate although strictly speaking, we expect deviations to exist. Once again, 
the above mentioned reasons for the relaxation rates come into play. Finally, we show that 
the extra 3D dispersion one expects for these bands in optimally doped YBa2Cu307 can 
resolve the lack of T dependence of the Knight shifts and O relaxation rate observed for this 
cuprate. 

With the same parameters, the Hall effect is calculated using standard Bloch-Boltzman 
theory and is shown to have the observed monotonic decreasing T dependence varying by 
~ 50% over the temperature range compared to the experimental 12 value of m 50%, with 
absolute values about 13 times larger than measured values. 

The linear resistivity is also obtained due to phonon relaxation across the bands. The 
sensitivity of and p to the doping follows naturally from the model. The qualitative 
features of the above experiments are straightforward from the considerations presented 
here. These features are very robust and do not depend sensitively at all on the choice 
of parameters used. This is particularly gratifying because these features are observed for 
a variety of different materials. Derivations of so many of the anomalous normal state 
properties of the cuprates lends confidence to the overall theory. 

Section II motivates and defines the IBP model. A Hamiltonian is written down and 
solved within the BCS mean field approximation. The orbital character of our bands is 
chosen. The optimal doping for superconductivity as the Fermi surface touching point is 
discussed. The lack of a unique phase for the attractive coupling matrix element Vkk' is 
discussed and the difference between Josephson tunneling in standard BCS-like intra-band 
pairing and IBP schemes is derived. IBP due to coupling through phonons is proposed. Up 
to this point the arguments and equations presented in Section II are quite general. There 
are several different channels through which IBP could occur. One particular channel is 
chosen and we restrict ourselves to singlet pairs only. There is nothing in the general theory 
that precludes IBP due to phonons with triplet pairs as far as we can see. 

The 4 key Josephson tunneling experiments 2,3 ' 13,14 are discussed. We show that with 
our choice of pairing channel, the two "d-wave" results of Wollman 2 et al and Tsuei 3 et al 
can be understood. The observed c-axis tunneling of Dynes 13 et al can also be understood as 
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due to BCS-like pairs formed away from the Fermi surface touching point. We are unable at 
present to decide whether the "s-wave" result of Chaudhari and Lim 14 supports or contradicts 
our model due to the complexity of the single particle tunneling matrix elements in this case. 

Section III presents quantitative calculations of the various normal state properties 
mentioned above and the final section is devoted to a summary of our results. 

II. INTER-BAND PAIRING MODEL 

A. The Model 

The T c of the cuprates is quite sensitive to the doping. Typically, there exists one 
optimal doping 15 for the highest T c and a rapid change in T c above and below this doping. 
A notable exception is YBCO 16 where one may argue that there are two optimal dopings, 
T c = 62K, 92K with an intermediate crossover regime. Regardless, whether we argue for one 
or two optimal dopings for YBCO, it is natural to associate the change in T c with doping to 
an attractive coupling that is strongly doping dependent. Experimentally, 17 many normal 
state properties are simultaneously strongly dependent on the doping. Thus, one is led to 
propose a bosonic excitation that couples to the electronic charge where either (or both) the 
bosons or their electronic coupling is strongly renormalized by the doping level. Rather than 
such dynamic schemes for the change in the attractive coupling with doping, we ask if there 
exists a simple kinematic scheme for the doping dependence of T c . In such a scheme, we 
hope to find a pairing strength that weakly depends on doping with the suppression of T c 
due primarily to a kinematic mechanism. IBP is such a mechanism. 

Consider two bands labeled U and L with dispersions ^\ ejj^ and further suppose 
there exists some attractive coupling leading to Cooper pairs of the form {kjj t, —ki I) 01 
(ki t, — ku l). In general, ej^ ^ e£\ It is only when ej^ and are both close to the 
Fermi energy that there is any possibility of a lowering of the overall energy by the formation 
of a Cooper pair. In BCS superconductors, pairs are formed for k states satisfying, 

|e* - cf| < fiwD, (1) 

where %uj]j is the Debye energy and e F is the Fermi energy. For IBP, the analogous criteria 
are 

\ef ] -e F \ <Ku, D , (2) 

1 4^ ~~ e F\ < Tu^d. (3) 

For BCS superconductors, there are always k states in the required energy range for pairing 
regardless of the doping. This is no longer the case for IBP because (2) and (3) must both 
be satisfied. It is only for special dopings that the U band and L band Fermi surfaces 
touch or are close to touching that any possibility exists for IBP. Hence, with IBP the 
doping sensitivity of T c can be understood quite simply as due to the change in the minimal 
separation in k space of the U and L band Fermi surfaces. No change in the strength of the 
pair attraction is required. In effect, we have replaced the requirement of a doping dependent 
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pairing interaction by a geometric argument based upon the band energy differences in k 
space. 

We will of course, take the mediator of the inter-band attraction to be phonons as in 
BCS. Because our pairs are comprised of electrons from two bands, the standard arguments 
for an upper bound of T c m 30K do not apply. This is certainly gratifying as it does not 
exclude IBP with phonons as a theory for high temperature superconductors. At present we 
have no estimate for the expected range of values for T c with our mechanism. 

There are some important differences between the pairing terms in the inter-band 
Hamiltonian and a BCS Hamiltonian. In order to highlight these differences we will briefly 
re-derive some of the well known aspects of the BCS Hamiltonian in order to both generalize 
to IBP and to establish our notation. 

In general, spin is conserved. Thus, there are two kinds of scatterings of BCS-like 
(k, —k) pairs: singlet to singlet, and triplet to triplet. Let the matrix elements for scattering 
a (k, —k) singlet (triplet) to a (A/, —k') singlet (triplet), 



-J= (0 fc 0_ fc ± 0- fc 0fc) (n t it) -+ -j= {<t>k><t>-k> ± <t>-k><t>k> ) -j= (n t tt) , (4) 

be V k \ (V k , k ). 4>k is a single particle wavefunction of momentum k. The first important point 
is that although <p k is uniquely defined only up to a phase <p k — > e l0 ^ with 9{—k) = —6(k), 
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V k , k is uniquely defined due to the time reversal invariance of a pair with itself. 

The Hamiltonian being Hermitian means V k , k = V kk l* and time reversal invariance of 
H implies V k , k = V kk '} . Thus, V k , k is always real. In second quantized form, 

-J=(0 fc 0_ fc ± 0_ fc fc )-^=(t; t It) = ^(4r a -fci ^ 4i a -fct)' ( 5 ) 

lie creat 
pair-pair interaction is, 

V k'k{l( a l'] a -k'i + a l'i a -k']-)( a -kl a k] + a-fcT°fc|) + 
k'k 

a \'^-k'] a -k] a k] + a k'l a -k'l a -kl a kl} + 

V k'k (g) (4'| a -fc'| - a l'i a -k'] )( a -kl a k] ~ a -k]a k [), (6) 
which can be expressed as, 

v k'k^\'^-k'i a -k[ a k] + j(Vk'k - V k '- k )al, a a^ k , a a_ ka a ka , (7) 
k'k k'k 



where a\ a is the creation operator for an electron in state <p k with spin a. The general 
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where, 



Vk'k = ^^k'k + V°k'-k + Vk'-k + V-k'-k) + ( v k'k + V-k'-k ~ V k'-k ~ V-k'k)]- ( 8 ) 

From the relations for V k °, k , we have V^k = V kk i m general, and V k 'k = Vk'-k for pure 
singlet pairing, and V k / k = —Vk'-k f° r pure triplet pairing. Finally, V k > k is always real and 
uniquely defined. For phonon coupling, V k / k is constant leading to pure singlet pairs. The 
symmetry of V^k determines the total spin of the Cooper pairs in BCS. This is not true for 
IBP. 

For IBP, let the matrix elements for the scatterings, 



-^{<l>Uk<t>L-k ± H-k<l>Uk)^={U T tt) -^{<t>Uk'<t>L-k' ± <i>L-w<i>Uw)-j=(U T tt), (9) 

be V k , k , V k , k . Here, 4>u,Lk are the single electron wavefunctions for each band with momentum 
k. In this case V k , k is no longer uniquely defined for a change in the definitions of <pu k and 

<pLk, 

<Puk -> e^W fcfc, (10) 
0Lfc^e^«0 Lfc , (11) 

where 9 UjL (-k) = -0 UjL (k) implies, 

y0,l _^ e -i[9u(V)-0 L (k!)] e i[eu(k)-9 L (k)] yO,l 
fc' 'k k^ k ' > / 

Hermiticity implies, = V^* and time reversal symmetry leads to V k , k = 

rather than V k , k = V kk l for BCS. This is due to inter-band Cooper pairs not being time 
reversal invariant with themselves. 

The general pair-pair inter-band pairing interaction is, 

V k\{~( a Uk'1 a L-k'l + a Uk'l a L~k'])( a L~kiauk] + O-L-k^Uki) + 

k'k 



a 



Uk'] " 



a L-k'] a L-k]CLUk] + a Uk' l a L-k' l a L-kiaUkl] + 



J2 V k'k ( (4/fc'T a I-fc'j _ a [/fc'| a L-/c'T)( a i-feI a ^T ~ «L-fcT a t/fcj)> ( 13 ) 

fc'fc V 7 

These equations make it clear that one cannot even speak of the symmetry of the pairing 
V^k without first specifying the phase convention on the single particle orbitals. Secondly, 
say we found a phase convention such that V k ,' k is pure s-wave. Unlike BCS phonon pairing, 



IBP does not preclude triplet pairing due to phonons. Experimentally, pure triplet pairs are 
not consistent with the observed Josephson tunneling between a BCS superconductor and 
the cuprates. There must be some singlet pairs for tunneling to occur. This does not exclude 
the possibility of some pairs being triplet paired, though. For the rest of this paper however, 
we will only consider the case of pure singlet pairs, i.e. V k \ k = 0. 
The IBP Hamiltonian becomes, 



e k a Uka a Uka + 4 a Lka a Lka 



ko 



+ 



n V k'k(a\jk'^ a L-k'i ~ a Uk'i a L~k'])( a L-kiauk] - a L - k }auki), (14) 
k'k z 

where we write V k i k = V k , k for simplicity. 

In the pairing approximation, we consider excitations of the form, 

Mk) = (u k + v k <J m a\_ k[ ) |0), po(k) = (1 - /f - fl L) ), (15) 

Mk)=a\j^\0), Pl (k) = fl U \l-fl L) ), (16) 

Mk) = al_ H \0), p 2 (k) = (1 - fl U) )fl L \ (17) 

Uk) = {-v k + u k a ] Uk] a{_ ki ) |0), p 3 (k) = & L \ (18) 

K| 2 + K| 2 = l, (19) 

with probabilities Pi(k) where f k U ^ and fj^ are the occupation numbers of U and L band 
particles of momenta k and —k respectively. There are four other excitations ■04,..., ipj 
obtained from the above by replacing o\j k ^a\_ k ^ with <^i_ k ^o\j k ^ in V'o an d ^3, a \jk} with 

a\_ k ^ in ip\ and with a\j k ^ in ip2- The ordering of the creation operators is chosen 

such that Cooper pairs consist only of singlet pairs. 

The coefficients u k and v k are chosen to minimize the free energy F = H — /iN — TS. 
The solution is, 

\u k \ = 1 - \v k \ = 2 I 1 + I , «fcVfc = (20) 

a = ^-/i, w * = i (ef } + ef } ) , (21) 
£ fc = ^ + |A fc | 2 , (22) 
^^fMV^ (23) 





4 '] + E k , 



(24) 



(U) 



1 



1 



(25) 



Ik 




Jk 



e^ +1 



with gap equation, 




(26) 



Conservation of the total number of particles iV leads to an equation for the chemical 
potential //, 



Unlike BCS, the quasiparticle excitation energies are different for the U and L particles 
and are given by and . These energies are the sum of half the difference in the U 
and L band energies and the term Ej, which is analogous to the BCS quasiparticle energy. 

The most important point of all is that, although the sum E^ + E^ = 2E k is always 
positive, the U or L excitation energy can be negative or have lower energy than the gap 
energy A. This new piece of physics will be considered in more detail below. 

Also, the BCS quasiparticle energy E^ is formed from a band that is the mean of 
the U and L bands. One can see why this is the case by noting that in the BCS ground 
state, pairs are either fully occupied or completely unoccupied. The individual band energies 

always appear summed together, ejjf 7 ^ + = 2u>k- For these states, the system "doesn't 
know" that the U and L pairing electrons have different energies. The three equations for 
the pair occupation amplitudes and Vk, incorporate the difference between the two band 
dispersions only through the gap This is to be expected by the same argument as above 
because Uk and Vk represent pair occupations. Similarly, the gap equation incorporates the 
U and L band differences through the quasiparticle occupation numbers, and f^ L \ 

Looking at the expressions for the quasiparticle excitation energies Ef ] and^ L) , one 

sees that the size of the difference in energies of the two bands ej^ — ej^ is what determines 
the size of the excitation. A negative energy implies that no Cooper pair is formed at T = 0. 
Instead, there is a single free electron in one band and no electron in the other. When 
e fc 7 ' ) ~~ e fc^ * s l ar g e > ^ is energetically unfavorable to occupy both k states or empty both k 
states. Thus no pair is formed. Near a Fermi surface touching point, ej^ ~ Therefore, 
inter-band pairs are always energetically favored. 

These equations are very similar to the equations for gapless superconductivity. 18 In 
gapless superconductivity, an spin electron has a slightly different dispersion than the \, 
spin electrons due to the presence of magnetic impurities. These t and \, spin bands are the 
analogues of the U and L bands in IBP. 



A r = £2M 2 (i-/f l -/f ) ) + (/f ) + /f ) ), 



(27) 



k 
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B. The Two Relevant Bands 

Perry and Tahir-Kheli 4 have calculated ab-initio the existence of two bands near the 
Fermi energy for Lai.ssSro.isCuC^. The results of that work are briefly summarized in this 
sub-section. 

La2Cu04 has two structural phases. A high temperature body-centered tetragonal 
crystal with point group and a low temperature orthorhombic lattice with point 
group. In the low temperature phase, the CuC>6 octahedra are tipped by 4.3° from their 
high temperature positions thereby reducing the symmetry of the crystal. The low temper- 
ature phase is the structure for superconductivity. In the tetragonal crystal, there are two 
reflection planes defined by the z-axis and the lines x = y and x = — y respectively. For the 
orthorhombic crystal, there is only one reflection plane defined by the z-axis and x = y. 

In either case, there will be a rigorous band crossing along k x = k y for our two bands. 
The tetragonal phase will also have a crossing along k x = —k y . The orthorhombic phase will 
come close to crossing along k x = —k y , but cannot cross. For IBP to occur, it is imperative 
that a crossing exist, for otherwise it is hard to see how IBP can overcome the band repulsion. 
For the remainder of this paper and in the following paper, we will take the LASCO crystal 
structure to be the high temperature tetragonal phase for simplicity. The small difference in 
structures can have a big effect on T c , but for most normal state properties, the difference 
will be small. 

In undoped La2Cu04, the O sites have a formal charge of —2, the Cu charge is +2 and 
is in a d 9 state. The point charge field on the Cu site due to the surrounding atoms make 
the d x 2_ y 2 orbital the most unstable leading to holes in this orbital in the undoped system. 
The next highest (unstable) orbital on the Cu is d z 2 due to the field of the apical oxygens. 
Undoped, each Cu d x 2_ y 2 has a single electron and d z 2 has 2 electrons. We propose that as 
the system is doped to its metallic/superconducting phase, holes start to appear in the Cu 
d z 2 orbital. At an arbitrary k point in the Brillouin zone, d z 2 can mix with d x 2_ y 2, O p a , 
and apical O p z orbitals. Thus, the two relevant bands at the Fermi energy are comprised 
primarily of these four orbitals. 

This assumption of the existence of d z 2 holes near the Fermi energy along with the 
assumption of inter-band pairing must be regarded as the two most important postulates of 
the IBP model for the cuprates. The former postulate is discussed in detail using ab-initio 
calculations on small clusters in the following paper 4 where we conclude that there are two 
bands near the Fermi energy with the character described above. 

We will call the two bands the Upper (U) and Lower (L) bands where eF^ > The 
lower band should be full or almost full in the undoped case and the upper band is half full 
with all of its holes predominantly of d x 2_ y 2 character. For LaSrCuO, the U and L bands 
will be almost completely 2D, while for optimally doped YBa2Cus07, the bands will have a 
measurable amount of 3D character. Underdoped YBa2Cu306.63 will have less 3D dispersion 
than optimally doped YBa2Cu307. Additionally, for YBCO there are potentially 4 relevant 
bands due to the 2 CuO planes per unit cell. Here, we restrict our attention to LaSrCuO 
and briefly discuss the differences to expect for YBCO where appropriate. 
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The unit cell of two formula units of La2Cu04 is tetragonal with lattice spacing a = 
4.0 A in the x, y direction and c = 12.0 A in the z-direction. The Brillouin zone of the 
primitive unit cell is given by — ir/a < k x < ir/a, —n/a < k y < ir/a, —2n/c < k z < 2n/c. 

To a first approximation, we may take the U and L bands for LaSrCuO to be purely 2D 
and add in the weak z-axis dispersion as a first order perturbation. This is done primarily 
for further computational simplicity and to convey the key aspects of the model. A more 
correct description will not qualitatively alter the behaviors we obtain for the NMR and Hall 
effect and more importantly, will not affect the general arguments for the various normal 
state properties. 

The relevant orbitals are d z 2, d x 2_ y 2 on the Cu, p a orbitals on the two planar O sites 
(a along the CuO bond direction) and p z orbitals for the two apical O sites above and below 
the Cu atom. Additional orbitals that appear are the Cu 4s, d xy , and O p n in the plane 
and the p x , p y orbitals on the apical O's. None of these additional orbitals will lead to a big 
change in the band structure but can affect the nmr. In particular, a small amount of Cu 
4s is required for an understanding of the sign of the Knight shift on Cu and some (a few 
percent) d xy is necessary for explaining the large anisotropy in the Cu spin relaxation rates 
for planar and z-axis magnetic fields. 

The tight-binding Hubbard model is (we use the electron picture), 

H = H ovh + #hop, (28) 
H or h = X! e d x 2_ y 2d x 2„ y 2 na d x 2_ y 2 na + ed z2 d) z 2 no d z 2 n(J + 

na 

e p CT (PxnaPxna + PynaPyna) + (pj/znaPUzna + PLznaPLzna), (29) 

na 

-f^hop = ^x 2 -y 2 ,a ^2_^2 ncr {Pxma + Pyma) ^tz 2 ,a ^z 2 na (Pxma + Pyraa) + 

{nm) {nm) 



^aa PxnaPyma ^ z tp z ,d z 2^2^z 2 na(P z Una PzLna) + ty z ,Vz ^^jP rJJnaPzLna + 

<nm>a na na 

^■taaa {PxnaPxma + PynaPyma) i ^p z ,a (p\ljna ~ PzLna)^Pxma ~ Pyma) + C.C., (30) 

{nm)a {nm) 
same a 
axis 

where d\ 2na creates a d z 2 orbital with spin a at site n. plna {Pyna) creates a p a electron on 

a planar O site along the x (y) axis and p\u na ^p\i na ) creates a p z electron on the apical O 
site above (below) the CuO plane, is the self (orbital) energy of the d z 2, etc. The t's are 
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the various hopping matrix elements. The ± sign in front of the hopping matrix elements 
represents the fact that the sign of the matrix element depends on the relative position of 
the two relevant orbitals. (nm) represents neighboring sites. The values of the parameters 
are shown in Table 1. The values in Table 1 for the hopping t's is for the antibonding 
combination of the two relevant orbitals. 

Except for the value of the difference 2 — ea 2 2 , there is nothing particularly surpris- 

z x —y 

ing or out of the ordinary with these parameters. €d 2 ~ e d 2 _ 2 being positive is what brings 
the d z 2 orbital up to the Fermi energy. This is the essential parameter for a two band de- 
scription like ours of superconductivity. These points are discussed in detail in the following 
paper 4 where our parameters are derived from ab-initio calculations on small clusters. As 
we stated in the introduction, the theory presented here is very robust and at a qualitative 
level does not depend sensitively on the values for these parameters at all. 

The point group of LaSrCuO is and d x 2_ y 2, d z 2 transform as B\ g and A\ g re- 
spectively. Under a reflection about the diagonals x = ±y, <Jdd x 2 ~y 2 = —dx 2 -y 2 an d 
<7dd z 2 — d z 2. Thus, for diagonal k vectors (k x = ±k y ) a single electron wavefunction must 
have d x 2_ y 2 character or d z 2 character. Along k x = ±k y , two bands may cross if one band 
has d x 2_ y 2 character and the other has d z 2 character. A plot of the top two bands (most 
unstable) is shown in figure 1 along the closed path (0, 0) — (it, it) — (it, 0) — (0, 0) in k space. 
For all k points not on the diagonal, there is no symmetry to keep d x 2_ y 2 and d z 2 from 

mixing. Thus, the two bands will repel leading to ej^ > e£\ 

The optimal doping for the highest T c is the doping that leads to a Fermi energy equal 
to the energy at the band crossing or touching point. It is at this doping and this doping 
only that there is favorable energetics for inter-band pair formation. Figure 2 shows the 
Fermi surface at various doping from underdoped to overdoped. 

We have adjusted the charge transfer of planar O p^ onto the La described in detail 
in the following paper 4 such that this doping matches the experimentally observed optimal 
doping for La2- a: Sr a: Cu04 at x — 0.15. Undoped, there are a total of 3 electrons in the U 
and L bands while for x — 0.15 there are 3 — 0.15 = 2.85 electrons in the two bands. 

The z-axis (normal to the CuO planes) dispersion is approximately calculated by ob- 
serving the dominant coupling in the z-direction is the p z orbital on the apical O above a 
CuO plane coupling to the p z orbital on the apical O below the next higher CuO plane. We 
take this value to be ^(0,0) = 0-3 eV. The coupling is added in via first order perturbation 
theory. Our choice for this matrix element is described below. 

The additional energy of a k state due to p z to p z coupling from layer to layer is, 

1 1 

e(k x , k y , k z ) = -2T( 0;0 )P cos -(k x a) cos -(k y a) cos k z c, (31) 

where P is the amount of apical O p z character in the k state (k x , k y ) derived from the 2D 
Hamiltonian (28). 

To incorporate the very small couplings through the p x and p y orbitals on the apical O 
above a CuO plane with the p z on the apical O below the next higher plane, an additional 
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energy, 



e(k x ,k y ,k z ) = -2T( 7r )P[cos-(A; :E a)sm-(A; 2/ a) + sin -(k x a) cos -(k y a)] cos k z c, (32) 
is added. Finally, we add, 

e(k x , k y , k z ) = -2T (?r ]7r) P sin -(k x a) sin ^(k y a) cos k z c, (33) 

to include the weak coupling of Cu d xy with apical p z . We take T(„. ) = 0.05 eV, 
T,,, : , 0.02. A'. 

The reason for naming the hopping matrix elements with the k space labels (0,0), 
(it, 0), (it, 7r) is that each one is multiplied by a combination of cos \k x a, cos \k y a, sin \k x a, 
sin ^A^a which is equal to 1 on the particular k space label and zero at the other two. We 
expect T (0>0) > T (7r>0) > T ( 

7T,7Tj • 

The primary effect of adding in the above z-axis couplings is to eliminate the logarithmic 
2D van-Hove saddle point singularity in the density of states of the U band by a broadened 
3D peak. The width of this peak is responsible for the T dependence (or lack thereof in 
YBa2Cu307) of the Knight shifts in the normal state. This is shown further along in the 
calculation of the nmr. 

Figure 3 shows the density of states of the two bands with the Fermi energy at the 
optimal doping indicated. The large peak in the U band density of states just above the 
optimal doping Fermi energy ep = 0.0 eV is due to the saddle point singularity at (ir, 0) and 
(0, 7r). The closeness of this peak to ep is a robust feature of this model and is not sensitive 
to the choice of parameters. 

Three other features of interest are shown in figures 4(a — c). They show the average 
amount of d x 2_ y 2, d z 2, and p a character of the two bands at different energies. Near ep, the 
two bands are predominantly d z 2. Both of these features are once again robust. 

C. The Pairing Term 

There are several different choices for the detailed orbital coupling that can lead to 
IBP. Lacking a rigorous microscopic proof that IBP exists, we make a particular choice for 
the pairing term and compare it to experiment. We take the pairing term in equation (14) to 
be mediated by an attractive coupling of the form shown in figure 5. In this figure, d x 2_ y 2] i 
and d z 2^ are defined as, 

<k»-y»k(r) = -^E eifc " fld *»-y*(r-R), ( 34 ) 



4 2k (r) = ^£ e ^M r -R)> (35) 

where R is the position of the Cu atom. Eventually, we will take the origin to be a center 
of inversion. Figure 5 depicts a d x 2_ y 2 Bloch orbital of momentum —k emitting a virtual 
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phonon and scattering to a d z 2 orbital of momentum —k' and a d z i orbital of momentum k 
absorbing the phonon and scattering to a d x 2_ y 2 orbital of momentum k' . (34) and (35) define 
the phase convention of the d x i_ y i and d z 2 Bloch functions and the matrix element Uk'k- 
The coupling Vk'k m (14) is calculated by projecting onto the above d x 2_ y 2, d z 2 scattering. 
Physically, what is accomplished is quite straightforward. The attractive coupling between 
two single electron states k, —k in the same band or two different bands is due primarily 
to the attractive coupling of the d x 2_ y 2 orbital component of one electron with the d z 2 
component of the other electron. One reason that such a coupling may be dominant is due 
to the d x 2_ y 2 k states localizing the electron charge in the plane, whereas the d z 2 —k state 
localizes the charge out of the CuO plane thereby reducing the Coulomb repulsion. 

An attractive coupling mediated by projection onto the diagram in figure 5 is also 
possible for (k t, — k I) pairs in the same band leading to a traditional BCS pairing term. 
Thus, one is faced with the possibility that BCS-like intra-band pairing may win out over 
IBP. This possibility is not realized for k states near the band crossing point by the same 
kinematics that allows a band crossing to occur along the diagonals (k x = ifcy). In order 
to have a large BCS-like pairing for a Cooper pair (k t, —k \) in the same band with a k 
state near the crossing point one requires a substantial amount of both d x 2_ y 2 and d z 2 in 
the k state. Near the band crossing point, the approximate diagonal reflection symmetry 
precludes that leading to a suppression of BCS-like pairing by the above pairing mechanism. 
Also, because the U and L k states are almost degenerate here, the system may further lower 
its energy by forming a pure d x 2_ y 2 and a pure d z 2 state in order to maximize the pairing 
coupling. At k points away from the diagonals, BCS-like pairing can occur and for arbitrary 
dopings, too. 

Including BCS-like intra-band pairing by projection of figure 5 as we did for IBP is 
straightforward. Rather than diagonalizing a 2x2 matrix as we did in deriving (15) — (27), 
we must diagonalize a 4x4 matrix. This is not difficult to do and leads to more complicated 
quasiparticle energies than (23) and (24). To properly include BCS-like pairs, we should 
not restrict ourselves to the single diagram in figure 5. There are other combinations of 
inter-band to BCS-like pair scatterings that may become important. 

Physically, what occurs is near the band crossing point on the diagonals, inter-band 
pairs are formed because one band is almost pure d x 2_ y 2 (with some O character) and the 
other band is almost pure d z i (again with O character). Away from the Fermi surface 
touching point, only BCS-like pairs are kinematically allowed leading to standard BCS pairs 
for the k states. Thus, one would expect to observe a gap in the photoemission for k 
vectors along the k x = and k y = directions. This is observed 19 for angle resolved 
photoemission (ARPES) on Bismuth 2212 although the observed Fermi surface crossing 
along TM ((0,0) — (n, 0)) has been questioned by Campuzano 20 et al as due to a k vector 
shifted by the superlattice translation (0.2l7r, 0.2l7r). 

We also expect a crossover region in the Brillouin zone where inter-band pairs change 
over to BCS-like pairs. Because IBP is only strong at the crossing point along k x = k y , the 
crossover region should occur close to the diagonals. For BCS-like pairs (k, —k), one would 
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see a gap A in ARPES at this k vector. On the diagonal at the touching point, one would 
again see a gap from (23) and (24). Slightly away from the diagonal while there are still 
inter-band pairs, (23) and (24) are still valid and at least one of these excitations has energy 
less than the gap. Thus, ARPES should see the "smallest gap" slightly off the diagonals 
k x = ±k y rather than on the diagonals. We have not yet calculated the size of this effect to 
determine whether present experimental resolution can see this effect. 

Owing to the separation of the pairs along the diagonal into one band with no d x 2_ y 2 
character and the other with no d z i character, pairs are most strongly formed here than at 
any other k point. This suggests that the dominant contribution to the Josephson tunneling 
current will come from IBP along the diagonals and we assume this is true unless such a 
current is rigorously zero. One can foresee that such a symmetry argument will be applied to 
inter-band tunneling in the c-axis direction making the contribution zero, thereby obtaining 
the "s-wave" result of Dynes 13 et al. 

Let, 

<p Uk = A Uk d x 2_ y 2 k + B Uk d z 2 k + other terms, (36) 

<\>Lk = A Lk d x 2_ y 2 k + B Lk d z 2 k + other terms, (37) 

where (f)jjk an d 4>Lk are the band wavefunctions and A^, B k are the projections onto the 
Bloch functions in (34), (35). Then, 

a Uk = A Ukd ] x 2_ y 2 k + B Uk d\ 2k + ■■■, 
4fc = A Lk dl 2 _ y2k + B Lk d\ 2k + • • • , 

and 

d x 2_ y 2 k = A* Uk a\j k + A* Lk a Lk H , 

d Uk = B Uk4jk + B Lk^Lk H • 

The pairing Hamiltonian term is, 

U k ' k (d x2 _ y2k ,^d z2 _ kli d x 2_ y 2_ ki d z 2 k ^ + d x2 _ y2k>) d z2 _ k ^d x 2_ y 2_ k ^d z 2 ki ). (42) 
Projecting onto (14), 

V k ' k = (A Uk >B L _ k ,)*(A Uk B L _ k )U k >- k + (B Uk < A L - kl )* (B Uk A L - k )U- k i k + 

(A Uk iB L _ k >)*(B Uk A L _ k )U k > k + (B Uk iA L ^ k /)*(A Uk B L _ k )U- k i- k . (43) 
Under the phase change in (10) and (11), 

Au, Lk -+ j QuAk) A v , Lk , (44) 
14 



(38) 
(39) 



(40) 
(41) 



Bu, Lk ^ e ie ^W B UjLk , 



(45) 



and thus V k 'k transforms as (12). 

Now every cuprate space group includes the inversion operator. We define the origin 
from which R in (34) and (35) is defined as a point of inversion. With this choice it is easy 
to see that the coefficients P^u,Lki ^>u,Lk must always satisfy A^B/j = real. 

In general, the electron-phonon matrix element that appears at each vertex in the 
Feynman diagram in figure 5 is the sum of terms of the form, 

i(e\ q ■ <l)V\ q (c q x + cl qX ) J e iq ' r d* x i_ y 2 k ,{r)d z 2 k {r)dT, (46) 

where A is the phonon polarization and q is its momentum. c q \ is the destruction operator 
and V\ q is the Fourier transform of the phonon potential and is the polarization vector. 
Inversion symmetry guarantees V\- q = V\ q , V\ q is always real and also that the integral in 
(46) over the electron states is real. 

Using (34) and (35), the integral in (46) is, 



/ d* x2 _ y2k ,(r)dMr)dr = / e*' £ e~ ik '- R ' e ik ' R d x 2_ y2 (r - R')d z 2{r ~ R). (47) 

RR' 

Taking the largest contribution to be when R' = R, 

J e^ r d x ._ y2k ,(r)dMr)dr = J e*' d x *_ y2 (r)d z z(r)dT. (48) 

Plugging this back into (46), we see that the product of the two electron-phonon matrix 
elements due to the two vertices in figure 5 is always positive with the phase convention 
defined in (34), (35) where the origin of R is an inversion center. With BCS-like pairing, 
time reversal symmetry of the Cooper pairs is sufficient to guarantee the product of the two 
vertices is always mod-squared and thus positive. 

Therefore, with our definitions of the d x 2_ y 2, d z i Bloch functions, we may take the 
total coupling U k > k from figure 5 to be s-like, 

Uk , k = I -V, \ef ' L) -n\< hw D and |e^' i) - n\ < hu D , ^ 
{ 0, otherwise, 

leading to, 

Vk'k = (-V)(Auk'BL~k' + Buk' A-L-k'Y (AukB L-k + B Vk A L _ k ). (50) 

The reason we have chosen the above range of k', k values for non-zero pairing is seen 
by considering the relevant Feynman diagram in figure 5. The value is, 
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{A Uk ,B L _ k i + B Uk iA L _ k i)*(A Uk B L _ k + B Uk A L _ k )M 2 , (51) 

where M 2 is a real positive number. If the constraints in (49) are satisfied, then the term in 
brackets is negative. 

D. Inter-Band Josephson Tunneling 

In BCS theory, a Cooper pair is of the form (f) k ^(f)_ k \ and transforms into itself under 
the operation of time reversal. For pairing across two distinct bands as we propose, the pair 
4>Uk} < t>L-ki transforms into (f)Lkl<f>U-kl and not into itself under time reversal. Of course, 
the full Hamiltonian remains time reversal invariant. This key difference alters the standard 
Josephson tunneling in a subtle yet dramatic way that leads to a new interpretation of the 
macroscopic phase of superconductors. 

Consider first the case of Josephson tunneling of a BCS-like Cooper pair on one side 
of a junction to a BCS-like pair on the other side of the junction. We will neglect all factors 
in the expression for the supercurrent J contributing only to the magnitude and not to the 
phase of J. The phase of the pair tunneling matrix element for the transfer of a (k t, —k \) 
pair to a (p t, — p I) pair is contained in the product of two single-electron tunneling matrix 
elements, T kp and T_ k _ p and the product of the two gap functions, A k and A' p on either 
side of the junction. This leads to supercurrent, 18 

J oc T kp T_ k _pA k A* = \T kp \ 2 A k A*. (52) 

In this expression, T kp is the matrix element for the transfer of a single electron of momentum 
k to an electron of momentum p, T_ k _ p is the corresponding matrix element for transferring 
the —k electron to — p and A k , A p are the gap functions on the two sides of the junction. By 
overall time reversal symmetry, T_ k _ p = T kp . Thus, the supercurrent is completely controlled 
by the phases of the superconducting gap functions on each side of the junction. The phases 
of the gap functions A, A' are determined by the symmetry of the pairing interactions for 
each superconductor. Thus, for all BCS-like pairing models, Josephson tunneling gives direct 
information of the symmetry of the gap. 

For IBP, the situation is dramatically different. Suppose we tunnel from an inter-band 
superconductor to a BCS superconductor. With Cooper pairing across bands, the matrix 
element for transferring a momentum k electron in the U band to a momentum p electron 

on the other side of the junction, is different from the matrix element for transferring 

a —k electron in the L band to a — p electron on the other side of the junction, T^ k . 

Although, = T^p* and = T^* by overall time reversal symmetry, the phase 

part of the pair transfer matrix element is of the form, 
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J*T^T^_ p A k A';. (53) 

In the case of inter-band pairing, the symmetry of the pairing interactions and the orbital 
character of the band wavefunctions both contribute to the overall tunneling phase. Under a 
redefinition of the single band wavefunctions (10), (11), 

^l^ierf)^^ (54) 

and from the gap equation (26), 

A fc e -i[&u(k)-e L (k)} Afcj ^ 

leading to no change in J. Note also that the above result is independent of the phase 
convention of the band orbitals on the BCS side of the junction because if an electron 
of momentum p is multiplied by a phase factor e % ®, then the phase of the —p electron is 
multiplied by e~~ %e . 

Remember that as we argued in the previous section, we take the dominant supercur- 
rent to be due to the inter-band Cooper pairs near the diagonal crossing points unless the 
contribution is rigorously zero. 

Let us apply our new relation (53) for the supercurrent to the four key Josephson 
tunneling experiments on YBCO. The first one by Wollman 2 et al on a YBCO-Pb corner 
junction is the simplest. We have shown that the phase of the supercurrent is independent 
of the choice of single particle wavefunctions. In this experiment, tunneling occurs along the 
two perpendicular Cu-0 bond directions in the CuO planes of YBCO (x and y axes) which 
are connected by a Pb wire. Pb is a BCS s-wave superconductor. A phase difference of tt 
implies a d-wave gap, whereas no phase difference implies an s-wave gap. A d-wave result is 
obtained. 

Let T^p'^ (x) be the matrix element for tunneling a k electron in the U or L band to a 
p electron in Pb along the x-axis. Choose the phase convention on the k states such that a 
90° (C4) rotation of the wavefunction for k is equal to the wavefunction for momentum C^k. 

Then, we must have T^F^ 4p (y) = T^ L \x). From (36) and (37), our phase convention 
gives, 

Au, L (C±k) = -A U:L (k), (56) 
Bu, L (C 4 k) = +B UjL (k). (57) 
Using the gap equation and the expression (50) for V k > k , 

A fc = A(A Uk B L _ k + B Uk A L _ k )*. (58) 

Hence, 



^C 4 k 



-A fc , (59) 



17 



leading to the result, 

J(C 4 k -> C 4 p) = -J(k^p). (60) 

Therefore, our IBP model also gives an observed "d-wave" gap. 

For the tri-crystal experiment of Tsuei 3 et al, we use the fact that for k near the 
Fermi surface touching point, (pjjk has almost no d z 2 character and 01^ has almost no d x 2_ y 2 
character or vice versa. For such a k vector, Bjjj- k, 0, A^ m (or Ajj^ ?s 0, m 0) and we 
can take the phase convention Ajjk, B^ real with > 0, B^ > (or Bjj^ > 0, A^ > 0). 
With this convention, is real and > 0. 

Suppose a d x 2_ y 2 orbital on one side of the grain boundary tunnels predominantly into 
a d x 2_ y 2 state on the other side of the boundary and thus d z 2 on one side goes predominantly 
to d z 2 on the other side. Then, the matrix element for d x 2_ y 2 — > d x 2_ y 2 is proportional to 
cos 28 cos 20 where 8 and are the orientations of the x and y axes in the CuO planes with 
respect to the grain boundary. Similarly, the d z 2 — > d z 2 matrix element has no orientation 
dependence and is thus constant. 

Suppose instead, that the dominant tunneling is d x 2_ y 2 — > d z 2 and d z 2 — > d x 2_ y 2. In 
this case, one matrix element is proportional to cos 28 and the other is proportional to cos 20 
leading the pair tunneling product cos 28 cos 20 as before. 

Taking the product of the various factors in (53), we see that the phase controlling the 
current is the same as to be expected from a d x 2_ y 2 gap as is observed. 

For the c-axis YBCO-Pb tunneling of Dynes 13 et al, an analysis similar to the YBCO 
corner junction case shows that the inter-band pairs do not contribute to the current. We 
still expect to see current due to the BCS-like pairs in this case. 

For the hexagonal YBCO tunneling of Lim 14 et al, the situation is different from the 
tri-crystal YBCO in one very important manner. Here a hexagonal MgO layer is placed 
on the LaA103 substrate and then YBCO is grown onto the sample. The hexagonal MgO 
causes YBCO to grow above it with its planar x, y axes rotated 45° from the angle of the 
YBCO grown directly over the LaA103 substrate. Thus, we expect a misalignment of the 
CuO planes across the six junctions complicating the matrix elements for single particle 
tunneling. Unfortunately, we have been unable to find a convincing argument telling us 
whether inter-band or BCS-like pairs dominate the single electron matrix elements in this 
case. 

III. NORMAL STATE NMR AND TRANSPORT 

A. NMR 

Armed with our parameters for the relevant bands and the criteria derived from IBP 
that the optimal doping is when the two Fermi surfaces touch, it is very easy to compute 
the normal state nmr and transport properties using standard expressions for the nmr and 
the Bloch-Boltzmann equation for the transport. It is quite satisfying that not only can we 
understand qualitatively the numerous anomalous features of the nmr and transport as due 
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to the character of the bands at this very special doping, but quantitatively the numbers are 
respectable. 

The key normal state nmr features 5-11 a theory for the cuprates must explain are: 1.) 
the difference in the relaxation rate curves with temperature T for the Cu and O nuclei in 
the plane, 2.) the similar Knight shifts (KS) at the two sites, 3.) the similarity of the O 
relaxation rate over T and the KS, 4.) the lack of T dependence of the KS for optimally 
doped YBCO7 and a monotonic increasing T dependence of the KS for optimally doped 
LaSrCuO and underdoped YBC06.63, 5.) the strictly monotonic decreasing Cu relaxation 
rate over T for optimally doped YBCO7 and the initially increasing and then decreasing Cu 
relaxation rate over T for underdoped YBC06.63 and optimally doped LaSrCuO, and 6.) 
the large anisotropy of the Cu relaxation rate for magnetic fields in the plane and along the 
z-axis. 

Figure 2 shows the two band Fermi surfaces at optimal doping. The electron-like band 
centered at k — (0,0) is the U(pper) band and the hole-like surface centered around (jr, tt) 
is the L(ower) band. Both bands are occupied at k — (0, 0) and unoccupied at (tt, tt). 

Figure 3 shows the density of states for the two bands. The first thing to notice is that 
at the Fermi energy ep, the L band has the lower density of states. We expect that this 
band predominantly carries the current and being hole-like will give a Hall coefficient with 
the correct sign. The large peak in the density of states of the U band at an energy a little 
larger 0.06 eV) than ep is due to the saddle point singularity at (tt, 0), (0, tt) for the U 
band. The width of the peak is due primarily to the overall strength of the z-axis couplings. 
This width is very sensitive to the details of the structure. This sensitivity controls the T 
dependence of the KS. The other aspects of the density of states are robust. 

Considering first the Cu spin relaxation, the two relevant orbitals are d x 2_ y 2 and d z i. 
Because there is no orbital relaxation between these orbitals regardless of the magnetic field 
direction, the relevant relaxation is dipole-dipole. We will neglect here the core polarization 
relaxation which we expect will not change the qualitative conclusions. The relaxation rate 
for a z-axis magnetic field is given by, 

™W X = 2(^)(7e7^) 2 / de/(e)(l - f{e))(^)\w^{e) + W& b ( e )], (61) 

^dip(e) = (^)[6Wd^(e)^ (62) 

W£ b (e) = 0, (63) 

where N{d x 2_ y 2){e) and N(d z 2)(e) are the total bare density of states at energy e and /(e) 
is the Fermi-Dirac function. 

The relaxation rate for a planar field Q3 W xy is identical to es W z above with W^ ip and 
W* Th replaced with, 
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(64) 




(65) 



Figures 4(a,b) show the amount of orbital character of d x 2_ y 2 and d z 2 in the U and L 
bands. The total bare density of states for d z i say, is the amount of orbital character of d z 2 
in the U band times the density of states of the U band Njj(e) plus a corresponding term for 
the L band. Figures 6(a — c) are plots of the bare density of states for each orbital due to 
each band and the total bare density of states. In all calculations of the nmr and Hall effect, 
the temperature dependence of chemical potential has been taken into account, /i decreases 
about 5 x 10~ 3 eV from T = to T = 300K. Finally, plots of 63 W Z /T and 63 W xy /T are 
shown in figures 7 (a, b) where we have used the value 21 (1/r 3 ) = 6 a.u. The curves show the 
characteristic peak at a T value greater than T c and the order of magnitude of the rate is 
consistent with experiment. 7 With the above figures and expressions for the relaxation rates 
we can qualitatively see why the correct behavior is obtained. 

There are two different ways an electron may relax a Cu nucleus: one by intra-band (U 
electron to U electron or L electron to L electron), or two, by the inter-band processes U to L 
and L to U. Due to the larger density of states of the U band, intra-band L to L relaxation is 
small. This leaves the relaxation rate to be determined by the U to U intra-band scattering 
and U to L (L to U) inter-band scattering. The U to U scattering leads to an increase of 
W/T due to the sharp increase in the U band density at the saddle point singularity. On the 
other hand, the contribution from the inter-band term must decrease due to the closeness of 
the Fermi energy to the very top of the L band and therefore the vanishing of the density of 
states. The competition of these two terms gives the final result shown. Regarding the Cu 
KS, only the U to U contribution and L to L contribution can appear because the KS comes 
from diagonal elements of the electron-nuclear Hamiltonian. The U to U KS contribution 
dominates the L to L KS contribution due to the larger density of states of the U band as 
before leading to a monotonically increasing KS if the hyperfine couplings have the correct 
sign. The sign of the couplings will be discussed later on, but we can now see that the Cu 
KS and W/T will have different temperature dependencies as observed. 

It is pleasing to find that our calculated value for the relaxation rate has the correct 
order of magnitude (experimentally, 63 W Z /T ?a 20 s _1 K _1 at 100K versus our value of 7.2 
s _1 K _1 ). Our calculation has completely neglected the effects of the Cu 4s contact term 
and secondarily, the core polarization contribution. Also, our calculated percentage increase 
from 30K to the peak is about 3% versus an experimental increase of m 10 — 20% and the 
percentage decrease from the peak value to the 300K value is 15% versus the observed 50%. 
In spite of these differences, the qualitative behavior is correct and this is the most important 
aspect given the level of calculation used in deriving the Hubbard model parameters. 
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As we can see from figures 7(a,b), the relaxation anisotropy of m 3.4 is not accounted 
for with the present model. What is missing here is a small amount of d xy orbital character 
in our bands. d xy is the next most unstable Cu orbital by ligand field theory after d x 2_ y 2 and 
d z 2. Including some d xy character affects W xy dramatically because now orbital relaxation 
is permitted W^ b ^ 0, whereas W* rb remains zero. Including d xy makes, 



There are also additional terms due to d xy in the dipole-dipole terms (62), (64) but these 
are small and cannot account for the anisotropy. We neglect them here. 

The coefficient of 4 in front of is about 30 times larger than the first coefficient 
in (62) and more than two orders of magnitude greater than the other coefficients in W^ ip . 
Thus, d xy character on the order of a few percent will lead to a contribution to 63 W xy as large 
as the dipole term leading easily to an anisotropy factor commensurate with experiment. 

For the planar O sites, the relevant orbital is p a . The O 2s will be considered later. 
The lack of a second major O orbital near the Fermi surface and also the fact that at (jr, ir) 
there is no mixing of antibonding (most unstable) p a (symmetry B l5 ) with d z i (A lff ) are the 
differences between O and Cu that leads to the different spin relaxation rate temperature 
dependencies. Figure 4c shows the amount of p a character on an O site at various energies. 
Figure 6c shows the p a bare density of states due to the U and L bands and the total bare 
density of states. 

We can see that the contribution to the bare density of states due to the L band 
A r L,j 3(T (e) is small compared to the U band term Nu tPtT (e) leading to a small contribution to 
the relaxation rate due to the L band. This is due to the L band k states near the Fermi 
energy being close to At (n,7r), d z 2 couples to the neighboring a orbitals which must 

be in a bonding configuration. This bonding set of a orbitals is stabilized by the t aa term. 
Thus, to create the most unstable d z 2 configuration at (ir, ir), we cannot have a large amount 
of p a character leading to the small value for N^ jPcr (e) near the Fermi energy. 

There are three distinct directions for the magnetic field at the O site with different 
relaxation rates. These are the z-axis normal to the CuO planes, the a axis along the Cu- 
O bond direction, and the perpendicular axis _L normal to z and a. With only the p a to 
consider, the relaxation along z and _L are equal, 17 W Z = 17 W±_ 7^ 17 W G . 

The expressions for 17 W Z , 17 W a are of the same form as for Cu (61) with the crude 
approximation (1/r 3 ) m 3 a.u. determined by ab-initio calculations and 7 n the gyromagnetic 
ratio of the 17 O nucleus. With only a p a contributing, there is no orbital relaxation in any 
direction of the magnetic field, 



K^) = ^ dx2 _ y2 (e)N dxy (e), 



(66) 



W&b = o. 



(67) 
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W£ b = 0. 



(68) 
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W^=(-^)N Pa (e)N Pa (e), (69) 

w£ p = wi = QwE p , (70) 

A plot of the 17 W a relaxation is shown in figure 8. The others are simply 2.5 times larger. 
This relaxation is monotonic increasing due to the bare density of p a due to the L band 
being small. 

There remains one orbital that can make a large contribution to the relaxation due to 
its large hyperfine coupling to the nucleus that we have not yet considered. That is the O 
2s orbital. A very small amount of 2s can have a large effect on 17 W without making any 
change in our calculated band structure. At (jr, ir), the L band has no d x 2_ y i character and 
by symmetry no 2s character either. Thus, similar to Ni p ^(e), there is almost no contribution 
to the O 2s bare density A2 S (e) due to the L band leading to the same monotonic increasing 
behavior as in figure 8. O 2s will make a large change to the absolute value of the relaxation 
rate. Without estimating the size of the 2s contribution, we cannot compare our answer in 
figure 8 to experimental data on LASCO. Without O 2s, the values in figure 8 are about an 
order of magnitude too small. 7 

The KS on the O sites is due to the O 2s Fermi contact interaction, the p a orbital 
dipole coupling and core polarization. In figure 9, we simply plot the KS, 

K a = -2K Z = -2K ± = Q(^y B J (-%)N p Ae)de, (71) 

due to the p a dipole-dipole term. Since the 2s and core polarization shifts are isotropic, our 
curve may be compared to the axial KS, = 1{K a — K±)/3. The monotonic increase in 
the shift is due to the sharp increase in the U band density of states above ep due to the 
saddle point singularity at (ir,0), (0, 7r). 

The KS involves only one factor of N Pa (e) while the relaxation contains N Pa (e) 2 . Thus, 
there is no way 17 W/T oc K a is strictly possible, although from the figures, the agreement 
is not far off. 

For the Cu KS, we need to introduce the effects of the Cu 4s orbital and the two p z 
orbitals above and below the Cu on the apical oxygens. Experimentally, the Cu KS for a 
z-axis field is T independent while for planar fields the shift is monotonic increasing with 
increasing T. As we discussed, the dominant contribution to the shift will come from the 
bare density of states for the orbitals from the U bands. If we can somehow show that the 
hyperfine couplings are positive for planar fields and zero for the z-axis field, then we should 
expect a curve similar to figure 9 for Cu KS. 

The Cu KS along the z-axis and planar directions is the sum of, 

63 K a = K^(d x ,_ y2 ) + K$»(d*) + iC ntact (4s) + + K^(d z2 ,As), (72) 
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where a is the field direction z or planar, K dip (d x 2_ y 2), K dip (d z 2) are the dipole shifts 
due to the d x i_ y i and d z 2 orbitals, ^ contact (4 s ) [ s the Fermi contact shift, K cp is the core 
polarization shift, and K dip (d z 2, 4s) is the dipole contribution due to the interference of d z 2 
and 4s. The expressions for the first four terms are, 

Kf p {d x 2_ y 2) = -2Kf p {d x 2_ y 2) = J (-^)N dx2y2 (e)de, (73) 

Kf p {d z 2) = -2K^(d z 2) = +*(^y B J (-^)N dz2 (e)de, (74) 

X contact (4s) = X contact (4s) = i^|^ s(0 )|V| / N 4s (e)de, (75) 

Kf = K cp = -{2a)Q^y 2 B J (-^[N dx2 y2 (e) + N dz2 (e)]de. (76) 

We take the value a = 0.33 (dimensionless) from Abragam and Bleaney. 21 ip4 S (0) is the value 
of the 4s orbital at the nucleus. 

The d z 2, 4s interference term is evaluated by taking the mean value of a band wave- 
function at the Fermi surface with the dipole Hamiltonian. Let, 

<Pk = A k d x 2_ y 2 k + B k d z 2 k + C k ij 4sk + . . • , (77) 
where d x 2_ y 2 k , d z 2 k are defined in (34), (35) and ip4 sk is defined similarly. Then, 

(0fc|#di P |0fc) = {-^B){lrfi)[I x S x + IyS y - 2I Z S Z ]- 

+ term involving (A* k C k + A k C£), (78) 




where I is the nuclear spin and R d 2 (r), i?4 S (r) are the radial parts of the d z 2 and 4s orbitals 
respectively with normalizations, 

J o + °° r 2 drR d Jr) 2 = ^ + °° r 2 drR 4s (r) 2 = 1. (80) 

Averaging over the Fermi surface, the interference term A k C k + A k C k due to d x 2_ y i 
and 4s becomes zero due to the different symmetries of d x 2_ y 2 and 4s (B l5 and A l5 ). The 
first two terms in (78) lead to K dvp {d x 2_ y 2) + K dip (d z 2) and the third term gives, 
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K^(d zh 4s) = (-^)(^) z2 /b J (-^)(B* k C k )N(e)de, (81) 

Kj$(d*,48) = -±K?>(d*,48), (82) 

where < B\Ck > is the mean value of B k Ck over the Fermi surface and N(e) is the total 
density of states of the band. Depending on the sign of < B^Ck >, the shift due to the 
interference term can be positive enhancing the field due to d z 2, or negative decreasing the 
net magnetic field of d z 2. The bare d z i and 4s density of states from this band are (\Bi t \ 2 )N(e) 
and (\C k \ 2 )N(e). Since only one power of C k appears in (81), (B^C^) can be large although 
(|Cfc| 2 ) may be small. 

The bare d x 2_ y 2 density of states is smaller than the d z 2 density of states at the Fermi 
energy leading to a net positive shift from K(d x 2_ y 2) + K[d z 2) along the z-direction and a 
net negative shift along the plane. x contact (4s) is isotropic and always positive. 

At all k vectors on the Fermi surface, the d z 2 and p z s on the apical O form an an- 
tibonding combination leading to a larger amount of 4s character than one expects from a 
d x 2_ y 2 band. Thus, if the d z 2, 4s interference term leads to a negative < B^C^ > that is 
sufficiently large, the net effect of K(d x 2_ y 2) + K[d z 2) + K(d z 2,As) leads to a dipolar field 
which is negative in the z direction and positive along the plane. We shall assume that 
due to the charge donation of the apical O p z to the Cu 4s, this is the case. With this 
assumption, we can see how the net shift on the Cu due to a z-axis field can add to zero 
while simultaneously leading to positive shifts for planar fields. In (72), the contribution to 
K a from the contact and core polarization terms should be a net isotropic positive value. 
Without the interference shift (81), there is no way the z-axis shift can add to zero. Figure 
10 shows the T dependence of the z-axis shift of K(d x 2_ y 2) and —K[d z 2). These curves have 
approximately the same T dependence as the O shift and 17 W/T. 

The above considerations also show why for fully doped YBCO7, the T dependence of 
the shifts and 17 W/T are almost constant. If there is more dispersion in the z-direction of 
the d z 2 orbital, then the saddle point peak in the U band density of states will be broadened. 
If the plateau is sufficiently broad to extend all the way to the Fermi energy or the (ir, 0) 
and (0,7r) saddle point singularity is sufficiently far away from the Fermi energy, then the 
shifts and O relaxation will become T independent. For optimally doped YBCO7, there are 
no vacancies in the chains and one can expect more 3D dispersion from the d z 2 orbital. 

B. Hall Effect and Resistivity 

Since the density of states for the L band is smaller than for the U band in the vicinity 
of the Fermi energy at the crossing point, we take the L band as the primary carrier of 
current. The Hall coefficient Rjj is the ratio of the transverse conductivity a xy and the 
conductivity a squared. Using standard Boltzman theory, 22 



fm 2 e n \fl\ / df\ 



d d 

Vk «dk x ~ Vk »dk y 



v kx , (83) 



24 



™a(^i;(-^r, (84) 



^=('^^. (85) 



where m e is the electron mass, fig is the primitive unit cell volume Qq = 96 A, f2 is the total 
volume, g = — |e| < is the electron charge, and c is the velocity of light. Vk = ^k^hl^ 1 is 
the velocity. We have multiplied a xy and a by the appropriate factors of m e and VLq to make 
the expressions in (82) and (83) dimensionless and have neglected the scattering rate 1/r in 
these expressions because for Rjj, r does not appear. We may also define <J xy (e) and er(e) 
by replacing the Fermi-Dirac function —df/de^ by the delta function 5(e^ — e) in (82) and 
(83). Then, 

a xy = J a xy {e)(^-^jde, (86) 

ff =.M-Hh (87) 

Figure 11 is a plot of cr X y( e ) an d cr(e) and figure 12 shows the temperature dependence 
of the Hall coefficient at optimal doping. The absolute magnitude of Rjj is about ten times 
larger than observed values for LASCO, but the calculated percentage change of Rjj from 
the 100K to 300K is ?s 40% in good agreement with experiment 12 (ps 50%). 

These curves were calculated taking the small z-axis dispersion parameters 7(o,o)> ^V,o)> 
and TJv,,.) in (31) — (33) to be zero. This was done to simplify the computation. <r xy (e) 
changes very rapidly for energies higher than the Fermi energy at the band touching point, 
whereas the change in <j(e) is not as dramatic. 

The abrupt change in u xy {e) is due to the two band crossing point. Near the Fermi 
energy of the crossing point, the repulsion of the two bands for non-diagonal k vectors 
strongly affects the shape of the two Fermi surfaces leading to a large curvature for the L 
band. At energies slightly higher than the Fermi energy, the L band surface is not affected 
by the presence of the U band and the Fermi surface becomes "small" leading to the reduced 
curvature. The overall T dependence of Rjj is due to the combined T dependencies of both 
a and o xy rather than solely due to o xy as one would expect from a first glance at figure 11. 
We can see that the presence of the two-band crossing point and the association of optimal 
superconductivity with this point is the root cause of the strong anomalous temperature 
dependence of the Hall coefficient. 

There are several possible explanations for our calculated value for Rjj being more than 
ten times too large. The most obvious one is the neglect of next-nearest neighbor hopping 
terms in our Hubbard model. The second cause is due to our 2D approximation to the band 
structure. A more detailed description than ours of the z-axis dispersion is required. The 
third reason is the neglect of the contribution from the electron-like surface of the U band. 
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Finally, more refined ab-initio calculations than the ones performed in the following paper on 
the Cu06 complex embedded in the point charge field of LASCO will change the parameters 
used in our Hubbard model. 

The resistivity due to scattering with phonons should be linear at the optimal doping 
for two reasons: 1.) the Fermi surface for the L band is "small", and 2.) at optimal doping 
the electron current can strongly relax by scattering to a U electron state. At optimal doping, 
the U band Fermi surface touches the L band surface leading to nearby states in k space 
with very different currents. The T dependence of a complicates this scenario because a is 
proportional to the effective number of charge carriers in the band. It is unclear how the T 
dependence of the relaxation rate 1/r could cancel this dependence. If a full-blown 3D model 
was used instead of our 2D model with the third dimension included as a perturbation, we 
believe that the temperature dependence of a would become small while a xy would remain 
very T dependent. 

CONCLUSIONS 

We have presented a model for superconductivity of the cuprates based on the idea that 
Cooper pairs are formed from electrons between two distinct bands. This leads naturally to 
associating optimal doping with a Fermi surface touching point of the two bands. We have 
postulated the character of the two bands to arise from d x 2_ y 2, d z 2, O p a , and apical O p z 
orbitals. A Hubbard model for these bands is set-up and we calculate some normal state 
consequences of the model. With our model, many of the anomalous features of the normal 
state nmr, Hall effect, and resistivity are explained qualitatively and to varying degrees 
quantitatively. The primary reason for the anomalous normal state properties is due to the 
optimal doping being at the Fermi surface touching point. 

We show that inter-band pairing alters the standard interpretation of Josephson tun- 
neling. With inter-band pairs, the detailed nature of the single particle tunneling matrix 
elements plays a prominent role. We show that with this new piece of physics, three of the 
four Josephson tunneling experiments are explained by our model with a phonon mediated 
attractive coupling. 

The parameters in the Hubbard model used in this paper are derived in the following 
paper by Perry and Tahir-Kheli from calculations on a Cu06 cluster for L^-xSr^CuO^ We 
conclude there that contrary to band structure calculations where only one band with Cu 
d x 2_ y 2 and O p a character is found, two bands exist at the Fermi energy with the character 
described above. 
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Table 1. Parameters in Hubbard Model in (eV). 



=2^03" 

e dz2 " -2.092 

t p -6.122 

e„ -0.852 

tx 2 -y 2 ,a 1-347 

t z 2 a 0.514 

ta„ 0.368 

t aaa -0.041 

t Pz4z2 1.076 

t Pz ,a 0.078 

t Pz , Pz 0.493 
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FIG. 1. Dispersions of the U (solid line) and L 
(dashed line) bands. 
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FIG. 2. Fermi surfaces at various dopings. The 
short dashes is for optimal doping, the long dashes is 
underdoping, and the solid line is overdoping. 
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FIG. 3. Density of states of the U (solid line) and 
L (dashed line) bands near the optimal Fermi energy. 
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FIG. 4a. d x 2_ y 2 character of the U (solid line) and 
L (dashed line) bands. 
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FIG. 4b. d z 2 character of the U (solid line) and L 
(dashed line) bands. 
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FIG. 4c. p a character of the U (solid line) and L 
(dashed line) bands. 
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FIG. 5. Feynman diagram for the pairing term 
with matrix element Uy^. 
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FIG. 6b. Bare density of states for d z 2 for the U 
(solid line) and L (dashed line) bands and the total 
bare denisty of states (dotted line). 
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FIG. 6c. Bare density of states for p CT for the U 
(solid line) and L (dashed line) bands and the total 
bare denisty of states (dotted line). 
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planar field. The contribution of d xy is not included. 
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